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Applying the method of integral estimates to the analysis of three-wave processes 
we derive the sufficient criteria for the hard loss of stability of the charged plane 
surface of liquids with different physical properties. The influence of higher-order 
wave interactions on the instability dynamics is also discussed. 

As we know , the dispersion relation for electrocapillary waves on a charged 

liquid surface has the following form: 

^ = g k+ -k 3 -—k 2 (1) 

where k = |k| is the wave number, u is the frequency, g is the acceleration of gravity, a is the 
surface tension, p is the mass density, and P is an external control parameter, depending 
on liquid nature. So, P = E 2 for an ideal perfectly conducting liquid (liquid metal in 
applications), where E is the external electric field strength. For an ideal dielectric liquid 
the parameter P equals the expression E 2 (e — l) 2 /(e 2 + e), where e is the permittivity. For 
liquid helium and liquid hydrogen it holds P = E + 2 + EJ 2 , where E + and E- designate 
the electric field strength above and below the fluid surface, respectively. One can see from 
Eq. (1) that if the control parameter P exceeds the critical value, 

P c = 8n^/gap, 

then uj 2 < and, as a consequence, an aperiodic instability develops. Thus, the condition 
P > P c is the criterion for the surface instability with respect to infinitesimal perturbations 
of the surface shape and the velocity field. 

As shown in Refs. , nonlinear interactions between three standing; waves, which 

form a hexagonal structure, can lead to the hard loss of stability of a charged liquid surface. 
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Then, even if the value of the control parameter is subcritical, i.e., P < P c , a fairly large- 
amplitude perturbation can remove the system from equilibrium, resulting either in the 
formation of a perturbed stationary surface profile or in infinite growth of the amplitudes 
of surface perturbations in accordance with character of higher-order wave processes. An 
important task of physical significance is therefore to obtain criteria for instability of the 
plane surface with respect to perturbations with finite amplitudes, namely, to determine the 
initial conditions (surface configuration and distribution of the velocity field) which lead to 
the development of an instability. 

In the present Letter we construct and discuss such criteria for the case of a conducting 
liquid in an external electric field (all our results may be extended to other fluids) by using the 
method of majoring equations, which was applied previously to the nonlinear Schrodinger 
equation (see, for example, Refs. 0, 0]), the nonlinear Klein-Gordon equation 1(3, 11 1, 
different modifications of the Boussinesq equation and so on. This method allows us 
to derive a number of sufficient integral criteria for the loss of stability of the plane charged 
liquid surface, with most of them corresponding to the subcritical values of the control 
parameter P, when the surface is stable in the linear approximation and the instability 
onset is caused by nonlinear processes. 

So, let us consider the irrotational motion of a perfectly conducting ideal liquid of infinite 
depth with a free surface z = r](x,y,t) in an external electric field E directed along the z 
axis. The velocity potential $ and the potential of the electric field (p obey the Laplace 
equations 



with the conditions at infinity: 



V 2 $ = 0, vV = o 



$ — > 0, z — > -oo, 



cp — ► —Ez, z — > oo. 

The conditions on the free surface are 

<9$ (V$) 2 (Vip) 2 -E 2 a„ V ± r] 

dri <9$ _ _ , 
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and, since the surface of a conducting liquid is equipotential, 



<f — 0, z = rj. 

The functions r](x,y,t) and <f)(x,y,t) = $| 2=r? are canonically conjugate [13], 
equations of motion take the Hamiltonian form, 

dip 8H dr] _ 5H 
dt Si] ' dt Sip' 

where the Hamiltonian coincides with the total energy of the system: 

H = ^kin + #pot> 



so that the 



(2) 



kin 



(V$) ! 



d 3 r, 



z<rj 



H 



pot 



d 2 r. 



(3) 



(4) 



Z>f] 



It is possible to express H explicitly in terms of the canonical variables. Rewriting the 
Hamiltonian in the form of a surface integral with the help of Green's formulas and expanding 
the integrand in powers series of ip and 77, we obtain 



H kin = / \ ( f +kf- l iP - V ±V ■ T+VlT - V) rf 2 r, 
#pot = " / + V ^ ■ ^-VjTr 1 ??) rf 2 r 



(5) 



d 2 r. 



(6) 



Here is the two-dimensional integral operator with a difference kernel whose Fourier trans- 



form is equal to the absolute value of the wave vector (ke l ^ r 
operators T± are defined by the expressions 



Ikle 1 ^ 1- ), and the nonlinear 



- ± 



E 

n=0 



{±rj) n k n 



nl 



For values of the control parameter e = (P — P c )/P c close to threshold (|e| <C 1), as is clear 
from the dispersion relation (1), the surface perturbations with k ~ ko — ^/gp/a have the 
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maximum linear growth rates. Then the nonlinear dynamics of the surface perturbations 
can be effectively studied with the help of the amplitude equation approach. Assuming 
the surface-slope angles to be small, |Vl^|, we introduce the slowly varying amplitudes Aj, 
j = 1, 2, 3, by means of the substitutions 

2 3 

V(r, t) = — > Ajixj^j, t) exp(zkj-r) + c.c, 

1 1 ,\ 2 ^— \ dAj(xj,yj,t) 
w — at — ex P^' r ) + cx -' 

j=i 

where the wave vectors k, with |k,-| = k make angles of 2n/3 with each other, and the 
variables Xj, yj form the orthogonal coordinate systems with Xj clXIS directed along the 
wave vectors k,. This representation for r\ and ip corresponds to a hexagonal structure of 
the surface perturbations, which is preferred at the initial stages in the development of the 
instability. 

Substituting the relations for rj and ip into Eqs. (2), (5) and (6), in the leading order we 
get the following expression for the averaged Hamiltonian: 

H = / E (l^l* + l^f - ^l 2 ) - a 1 a 2 a 3 - a{a;a; d\ (7) 

J U=i 

where we convert to dimensionless variables, 

r^r/( v / 2A; ), H — > 8gH/ (9/uq), 

and introduce the linear differential operators 

did 2 

The equations for the complex amplitudes Aj corresponding to this Hamiltonian are 

q2 a . 4*4*4* 

^ = eA j + L^A j + ^^, 3 = W- (8) 

It should be noted that the multiplier p — (e — l)/(e + 1) and the multiplier s = (E + 2 — 
E_ 2 )/P c appear before the nonlinear terms in the right-hand sides of Eqs. (8) for dielectric 
liquid and for liquid helium, respectively (the equations reduce to the form (8) by scaling 
the amplitudes Aj — > Aj/p or Aj — > Aj/ s). 
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Let us find the sufficient criteria for unlimited growth of the amplitudes Aj in a finite 
time, i.e., the criteria of the blow-up for Eqs. (8), which describe the nonlinear interactions 
between three standing waves. Consider the time evolution of the following positive quantity: 

3 



X = J2X„ X J (t)= ! \A \ 2 d 2 r. 

7=1 J 



Differentiating X twice with respect to t and then making use of Eqs. (8), we get 



. 3=1 



2 

3 It 



LjAjl 2 + e|A,-| 2 ) + 3A X A 2 A 3 + 3A{A* 2 A; 



d 2 r. 



Excluding the cubic terms from the integrand with the help of the expression (7), we come 
to the relation 

3 „ 

X tt + 3H = -eX + J2 [ 5 l^'l' + l4Af] d2r - ( 9 ) 

3=1 J 

It follows from the integral Holder inequality for the functions \Aj\ and \Aj\ t that 

AX 3 j \A 3 \ 2 t d 2 r > X 2 . 
On the other hand, the algebraic Cauchy inequality yields 

As a consequence, we can estimate the second term in the right-hand side of Eq. (9): 

' - x t 2 



3=1 J 



AX 



In addition, taking into account that 



\LjAj\ 2 d 2 r > 0, 



we obtain the following second order differential inequality: 



X tt + 3H > ^y-eX. 



(10) 



Solving the inequality one may find the sufficient conditions under which the integral value 
X becomes infinite in a finite time. It should be noted that the similar majoring inequalities 



were derived in Refs. |8l E 



10 



11 



as a result of investigation of the blow-up in different 



well-known nonlinear partial differential equations. 
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Under the substitution, Y = X 1(/4 , inequality (10) takes the Newtonian form 




(11) 



where Y can be considered as a coordinate of some particle and the function P(Y) plays the 
role of the potential. Let the particle velocity Y t be negative (in this case X t > 0). Then, 
multiplying the inequality (11) by Y t and integrating it over time, we get 



i.e., the particle energy U increases with time. It is clear that if the condition JJ t = 
holds, which corresponds to the equality sign in the expression (11), and the particle does 
not encounter a potential barrier, then it reaches the origin, i.e., the point Y — 0, and, 
consequently, the positive-definite quantity X becomes infinite. The collapse takes place 

(a) if e < and H > 0, provided that Y(t ) < |ej s /(3)s and l2U(t ) < \e\*/(3H)*; 

(b) if e < and H > 0, provided that 12U(t ) > \e\*/(3H)*; 

(c) if e < and H < 0; 

(d) if e > 0, provided that U(t ) > 0, 

where t — t corresponds to the initial moment. The collapse time t c at which the amplitudes 
Aj go to infinity can be estimated from above as follows, 



Note that the condition Y t (0) < is optional for the cases (a) and (c). Since U t (t) < for 
Y t > 0, in these two cases the particle always reaches the point Y = after the reflection 
from the potential barrier. 

It is important that the conditions (a)-(d) may serve as the sufficient criteria for the 
instability of the plane surface of a conducting liquid in the near-critical electric field with 
respect to finite amplitude perturbations. This distinguishes our criteria from the simplest 
criterion for linear instability, P > P c , which corresponds to infinitesimal perturbations. 
Actually, the conditions (a), (b) and (c) are relative to the case of the subcritical electric 



Utit) > 0, 



U(t) =Y t 2 /2 + P(Y), 



Y(to) 
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field strength, P < P c , when the surface is stable in the linear approximation. It is obvious 
that we deal with hard excitation of the electohydrodinamic instability. 

Thus, we have shown that if the conditions (a)-(d) are valid, then the equations (8) 
describe the unlimited growth of the amplitudes Aj. However, the applicability of Eqs. (8), 
i.e., the possibility of limiting the treatment to three- wave processes, assumes that the 
perturbation amplitudes are small (\Aj\ are the order of magnitude of the control parameter 
e). The question arises as to whether the neglected higher-order nonlinearities stabilize the 
instability, or, on the contrary, higher-order wave processes promote an explosive growth of 
the amplitudes. Notice that both the experimental data ^| and the results of the numerical 
calculations [15] indicate that the higher-order nonlinearities have the destabilizing influence. 

Let us show that, in particular, the four- wave interaction does not saturate the explosive 
instability. Consider the simplest case when the perturbation amplitudes do not depend on 
the spatial variables x and y, and the surface configuration is given by 

3 



V (r, t) = + 2 ^J) + ^(5v/3 + 6)(f i£ + + 6£) + cc, 

3=1 



£j{r,t) = a i (t)exp(ik i r), 

where the nonlinear interactions between fundamental and combination harmonics, ko <-> 2fc 
and k <-> V3k , are taken into account. Substituting these expressions into (6), we obtain 
the fourth-order correction for the potential energy (4): 

3 



<ot = -^o£ / ^K-| 4 + i(25V3 + 13)K.|- 2 n 

3=1 J L 3=1 



\ a j 



1 2 



d 2 r. 



This functional is the negative-definite quantity, < 0. Consequently, we can assume 

that the potential energy Hp Q ^ decreases indefinitely as the perturbations grow. Then 
the kinetic energy #v m , which is the positive-definite quantity, increases infinitely (see 
Eqs. (2) and (3)). Thus, the higher-order nonlinearities do not retard the explosive growth 
of the amplitudes in the model (8), and the integral criteria (a)-(d) can be considered as 
the sufficient criteria for the unlimited growth of the perturbations of a conducting liquid 
surface in an applied electric field. 

As for a dielectric liquid in the near-critical electric field, we have 

3 « r / ^\ 3 



<ot = -^o£ / (V - ~) M 4 + \ (15 -12V3+ P 2 (37v/3- 2)) \a^\[ 

3=i J y ' 3=i 



a i\ 



d 2 r. 
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It can readily be seen that < for arbitrary aj if e > E2D holds, where E2D ~ 3.53 

(at this value of the permittivity the hard instability regime changes to the soft one in 
2D geometry ^^). Then the relations (a)-(d) represent the criteria for the blow-up- type 
dynamics of the surface perturbations. If 1 < e < e^d holds, where e^d ~ 2.05 (the 
functional with | | changes the sign precisely at this value of the 

permittivity e [6]), then > and, consequently, the four- wave processes can stabilize 

the instability, resulting in the appearance of the stationary hexagonal structures. In this 
case the conditions (a)-(c) are the criteria for the hard excitation of the stationary wave 
patterns on the free surface of an ideal dielectric liquid. 

For liquid helium (or hydrogen) with a charged surface it holds 

3 „ r / k\ 3 



<ot = -^o£ / (4, 2 -^K| 4 + i(l6v / 3-29 + S 2 (9V3 + 42))K.|- 2 n 
j=i J V \ / j=i 



\ a i\ 



d 2 r. 



One can find that < for \s\ > \s2d\, where sojj ~ 0.56 (at this value of the parameter 

s the functional changes its sing in 2D case U), and H® t > for |s| < | s 3 d|, where 

s 3D « 0.33. 
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